


Roots by Iteration 



Welcome to another install- 
ment of the Programmer's 
Toolbox. Lately, we have 
been looking at algorithms to find the 
root of a function: 



y = f(x) 



(1) 



Two months ago, we looked at 
Newton's method, which is the first- 
order method of choice if you have an 
analytical expression for the derivative 
of f(x). Newton's method is fast but 
unreliable unless we take steps to tame 
it by limiting its range of search. Last 
month, we looked at three methods we 
can use if the derivative is unavailable: 
the method of bisection, the secant 
method, and the method of false posi- 
tion. The last two methods converge at 
roughly the same rate as Newton's and 
m are called super-linear methods. 

We learned that each method has its 
pros and cons. The bisection method is 
like the tortoise, slow but sure. It is 
simple and converges at a guaranteed 
pace. If the function is continuous and 
has a root in the region of interest, this 
method will find it, and we can even 
predict how many steps it will take to 
do so. 

The secant method converges faster 
than competing methods, but like its 
philosophical equivalent, Newton's 
method, it can also oscillate or diverge. 
Special steps should be taken to tame 
the method if it is going to used in a 
general-purpose routine. The method 
of false position is much more stable 
and, as with the bisection method, we 
can guarantee convergence. Most of 
the time, convergence is as fast, or 
nearly so, as the secant method. But for 
some cases (and not just pathological 
cases either), the convergence can be 
much slower — even worse than the 
method of bisection. 

Throughout this study, I've warned 
you to be careful applying any of these 



A first-order 
method is 
on the idea of 
fitting a straight 
line between two 
points lying on the 
function f(x). 

algorithms blindly— to "know your 
function." Without care in implemen- 
tation, all the algorithms can lead to 
trouble. We need a method that's more 
trustworthy, a method that combines 
the stability and trustworthiness of the 
method of bisection with the fast con- 
vergence of a secant or Newton's 
method. This month, I'm going to do 
even better: I'll show you a second- 
order method that's not only stable and 
trustworthy, but converges with light- 
ning speed once we get near the root. 

METHOD X 

The algorithm I'm going to show you 
here has an interesting history. I first 
encountered it in the old IBM 
Scientific Subroutine Library, a library 
of Fortran functions, circa 1965. IBM 
described the algorithm as Mueller's 
method. It's based on the concept of 
inverse quadratic interpolation, which 
I'll explain more fully in a moment. 
Many years ago, however, I saw a 
description of Mueller's method and 
discovered that it uses direct quadratic 
interpolation — a critical difference, as 
you'll soon see. Numerical Recipes in 
C describes Brent's method, which 
also uses inverse quadratic interpola- 
tion and has many of the same charac- 



t - stics as the IBM algorithm. 1 But the 

used to decide whether to accept M 
approximation are all different. That 
leaves me with a method that works 
very well, but no one to attribute it to. 
If any of you can identify the inventor 
of the algorithm, please let me know. 
Until then, I'll just call it Method X. 
The derivation of the acceptance crite- 
rion is my own. 

A rose by any other name would 
smell as sweet, and Method X is sweet 
indeed. It shares the property of the 
method of bisection: If a root exists in 
the range specified, this algorithm will 
find it— guaranteed. But unlike the 
method of bisection, it converges 
extremely rapidly. If the secant method 
is super-linear, this one is super-qua- 
dratic and will find the root so fast it 
will make your head spin. Who could 
ask for anything more, except a name? 

THE BASICS 

A first-order method is based on the 
idea of fitting a straight line between 
two points lying on the function f(x). 
Not surprisingly, a second-order 
method is based on fitting a quadratic 
through three points. 

Consider the curve of Figure 1 . We 
can fit a quadratic through the three 
points shown to get the parabola 
described by the dashed curve. The 
general form of such a quadratic is: 



y = a + bx + cx 2 



(2) 



but by playing around with the coeffi- 
cients, we can eliminate the first-order 
term to get: 



y-y =A(x-x o y 



(3) 



This equation describes a parabola 
opening up (for positive A) or down 
(for negative A), with vertex at the 
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point (x , y ). By looking at the values 
of the function at x = x,, x 2 , and x 3 , we 
could conceivably determine the con- 
stants x , y , and A. The root is then 
given by setting y = to get: 



(4) 



This is Mueller's method. Now, look- 
ing at Equation 4, tell me: what's 
wrong with this picture? That's right, 
we have two roots (potentially com- 
plex), not just one. The function is dou- 
ble-valued when solved for x as a func- 
tion of y. How do we decide which of 
the two roots is the right one? And 
what if the roots are complex? In that 
case, we get no useful solution. 

While we can conceivably make 
tests on the results to be sure Equation 
4 leads to the right root (we'll end up 
performing similar tests, in any case), 
wouldn't it be nice if our algorithm 
could assure us one real root in all 
cases? 

This is where the inverse part of 
"inverse quadratic interpolation" 
comes in, and with it, the ingenuity of 
Method X (and Brent's method). The 
concept is very simple. Instead of a 
parabola opening along the y-axis, let's 
use one that opens along the x-axis; 
that is, it will be a quadratic in y rather 
than x. A function like this will always 
be single-valued, since each value of y 
gives us one and only one value of x. 
Since y can take on any value, we 
needn't worry about complex or imag- 
inary roots. An inverse quadratic fit is 
shown in Figure 2. 

When you're trying to fit a curve 
through known points, the ease of the 
process depends a lot on the form you 
choose for the function. A general 
form like Equation 2 may look simple 
and straightforward, but you'll find 
that when you start substituting in the 
values at the known points, you're in 
for some tedious and error-prone math. 
I've found a trick that often simplifies 
things quite a bit. Write the equation in 
a factored form: 

x-x, = A + B(y- yi )[\ + C{y-y 2 7l 

(5) 



You can see the advantage to this 
form when you evaluate it at the 
known points. At point Pi, we have: 



x=x, 

y=y i 

FIGURE 1 

Quadratic Jit. 



(6) 



so the equation simply gives us: 

0=A 

Thus, A must be zero. 
The next coefficient is equally easy 




FIGURE 2 

Inverse quadratic fit. 



y = f(x) 
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since, at point P2, the equation 
becomes: 

x 2 -Xi= B(y 2 - yi ) 
This gives: 



C = 



(7) 



As you can see, we've already deter- 
mined two of the three coefficients 
with very little work. We get the last 
coefficient by evaluating at P 3 : 

x 3 - *, = B(y 3 - y,)[l + C(y 3 - y 2 )] 
Solving for C, we get: 

Using Equation 7, this becomes: 

-x\){yi -y\) 



l + C(y 3 -y 2 ) = 



(x2 -xi)(y 3 -yi) 



(8) 



So far, we haven't placed any 
restrictions on the values of x b x 2 , and 
x 3 . But let's now assume that x 2 is 
halfway between x, and x 3 , so: 



x 2 — x [ — x 3 — x 2 — h 
Equation 7 becomes: 

Equation 8 becomes: 



(9) 



00) 



i + c(y 3 -y 2 ) = 



2(y 2 -yQ 
y 3 -yi 



Continuing with the solution for C 
gives us: 

V/3 72/ ^ _ yi 

= 2(y 2 -yi)-(y3 - yi) 
ys - yi 

_ -yi + 2y 2 - y 3 

ys - yi 

and thus: 



-y, + 2y 2 - y 3 

(ys -y 2 )(y3 - yi) 



(ID 



We now have all three coefficients in 
terms of the fitted points. Remember 
this "factored polynomial" trick — it's 
very useful, and works for polynomials 
of any order. As you can see, the com- 
putation of the coefficients gets pro- 
gressively harder, but it will still be 
easier than with any other form. 

Once we have the coefficients, we 
can find the desired root by setting 
y = in Equation 5: 



x = xi - By,(l - Cy 2 ) 



(12) 



At first glance, you might think we're 
home free at this point, since we need 
only solve Equation 12 for x. 
Unfortunately, a lot of snakes still lurk 
in the equations to bite the unwary. For 
starters, if any pair of y-values are 
equal, one of the terms in the denomi- 
nators of Equations 7 or 1 1 will go to 
zero, and our equations blow up. I was 
very careful not to show the inverse 
quadratic fit in Figure 1. That's 
because y 2 = y 3 in that figure and no 
(finite) parabola opening horizontally 
is possible. 

FIGURE 3 

A poor interpolant. 



By requiring that points P] and P 3 
straddle the root, we can be sure that: 

yi*y3 

So at least one term can't vanish. But 
we have no control over the value of 
y 2 , so we can't promise that one of the 
other two terms won't vanish. What do 
we do in this case? That's easy: We 
simply refuse to do the interpolation. 
In finding y 2 , we have, in effect, taken 
one step of the bisection method. If we 
find that we can't do the interpolation, 
all is not lost. We simply take another 
bisection step and repeat the process. 
Eventually, the method must give us a 
useable interpolant. 

Even if the terms in Equations 7 and 
1 1 don't blow up, we might still find 
that our interpolation formula doesn't 
help. Consider the curve of Figure 3. 
While our formula works, it predicts a 
root that's well outside the range 
X; .. x 3 . Since we know that the desired 
root is in this range, we'd be foolish to 
accept the interpolated value. 

ACCEPTANCE CRITERION 

We need a nice mathematical test to 
tell us when we can trust the root esti- 
mated by the inverse interpolation. The 
obvious approach is to simply look at 
the predicted value of x. If it's not 
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inside the range of our search, we 
should refuse to accept the step. That's 
what's done in Brent's method. 
Unfortunately, we can get into trouble 
even then. Look at Figure 4. In this 
case, we get a mathematically correct 
fit to the data, and the predicted root is 
in the range of interest. But the accura- 
cy of the fit is terrible; the predicted 
root, near P 3 , is far from the actual one 
near ? v We'd have been better off just 
accepting a bisection step. What's 



more, if we carry the method to the 
next step, we'll find that the next pre- 
dicted root will be out of range. 

Deriving a suitable acceptance crite- 
rion is tricky, and I spent a long time 
finding it. But the final criterion is 
quite simple, and easily evaluated in 
software. The derivation is also fasci- 
nating in its own right, so bear with me 
as we delve into some tricky math. The 
results will be worth it, I promise. 

We start by assuming that the root 



lies between points P, and P 3 , and nei- 
ther point is the root (if it is, we're fin- 
ished). This assumption means that y! 
and y 3 have different signs and so: 

viy 3 < 

After bisection to get y 2 , we can tell 
from its sign whether the root lies to 
its left or its right. This is the basis, of 
course, of the method of bisection. 
Naturally, we don't know which of 
the two cases we'll get, but I wish we 
did. Looking at Equation 10, we can 
see that the relationship between the 
three points is not symmetric. That is, 
points P! and P 3 enter into the equa- 
tion much differently, thanks to our 
factored form. To avoid having to 
deal with multiple cases, we would 
like very much always to know that 
the root lies between P| and P 2 . We 
can assure this situation by the simple 
expedient of relabeling the points, if 
necessary. 

If the root lies between P, and P 2 , yi 
and y 2 must also have opposite signs, 
and y 2 and y 3 must have similar ones. 
Thus: 

viv 2 <0 

v 2 v 3 >0 (14) 

Now we'll assert our acceptance cri- 
terion: the predicted root must also lie 
between x, and x 2 . This is equivalent to 
the assertion: 

(x 2 -x)(x- Xl )>0 (15) 
From Equation 12: 

x 2 -x = x 2 -[xi- B yi (l - Cy 2 )] 
= x 2 -x, +fiy,(l-Cy 2 ) 

or: 

x 2 -x = h + fiy, (1-Cy 2 ) (16) 

(Because of our relabeling, h may no 
longer be positive.) 

Substituting for B from Equation 7, 
we get: 
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FIGURE 4 

A poor fit. 




x 2 - x = h + 
1 + 



4 



(1-Cy,)(l-Cy 2 )>0 



(19) 



y2-y\ 



or: 



x 2 - x 



yi ~y\ +yiO-Cy 2 ) 
yz-y\ 



Next, we must substitute the value 
for C from Equation H. After some 
truly horrible algebra that I'll spare 
you, we find: 



Similarly: 



(18) 



1 - Cy, = 



And: 



1 - Cy 2 = 



Let: 



yijys - yi) - y\{>-2 -yi) 



-yi)(y3 -y\) 



(20) 



yifa -yi)-2y 2 (y 2 -yQ 



(ys -y2)(ya - yi) 



Then the inequality in Equation 15 
becomes: 



-ft 2 y,y 2 



(1-Cy,)(l-Cy 2 ) 

(y2-yi) 2 



>o 



Now, h 2 and the denominator are, of 
course, perfect squares, so their signs 
won't affect the outcome. Also, we 
know from Equation 13 that the sign of 
y^ 2 must be negative. This leaves us 
with the condition: 



M = y3(y3-y2)-yi(y2-yi) 
v = y3(y3-yi)-2y 2 (y 2 -y,) 

Then we have: 
1 - Cy, = 



(21) 



(22) 



1-Cv 2 = 



(ys -y2)(y3 - yi) 



(>3 -y2)(y3 - yi) 



(23) 



When we multiply these equations 
together, we again will get perfect 
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squares in the denominator, which 
won't affect the sign of the result. So 
our inequality becomes, finally: 



mv>0 



(24) 



We can identify four cases, depend- 
ing on the signs of u and v: 



• Case 1: u > 0, v > 

• Case 2: u > 0, v < 

• Case 3: u < 0, v > 

• Case 4: u < 0, v < 



OK 

Forbidden 
Forbidden 
OK 



PUSHING THE ENVELOPE 

Now, we must somehow relate these 
cases to the corresponding values of 
the function at points P,, P 2 , and P 3 . 
We can explore the boundaries 
between the cases by setting u and v, 
separately, equal to zero. For example, 
when u = 0, we have: 

B(B -»)- »fe -») 54 ( 25 ) 

This equation is quadratic in yj and y 3 , 
but only linear in y 2 , so we can solve 
for that parameter. We get: 

y 2 (yi + y3) ■ yi 2+ y3 2 



or: 



y 2 = 



y. 2 + y.3 2 
y\ +y3 



(26) 



At this point, let's introduce two 
new variables: 



' y, r y x 



(27) 



Because of our restrictions on thei 
ranges (see Equations 13 and 14), both 
P and y must be negative. 

Using these definitions, we can 
rewrite Equation 26 in the form: 



1 + y 



l + y 

Rearranging the terms gives: 



,2 _ 



fiy + 1 - = 



(28) 



(29) 



This is a quadratic in y, which we cai 
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solve using the quadratic formula. Let: 
a = \; b = -p; c = 1 - jS (30) 
Then the quadratic formula gives: 



becomes: 



_ /?±V^-4(l-/Q 
Y 2 



(31) 



This curve is plotted as the lighter 
curve in Figure 5. As you can see, it 
has two lobes, enclosing the two areas 
marked "Case 4." 

We can get the other boundary by 
setting v = 0. From Equation 22, this 
gives: 



y 3 (y3-yi)-2y 2 (y 2 -y,) = o 



(32) 



This time, it's y, that we can solve for: 



yi ~ y 3 -2v 2 



(33) 



1 



y 2 - 2fi 2 
~ y-ip 

A bit of rearranging gives us: 

y - 2/3 = y 2 - 2p 2 

or: 

2/? 2 - y 2 - 2/3 + y = 



(34) 



This equation is a little different from 
the previous one because it's a qua- 
dratic in both p and y. In fact, it repre- 
sents a hyperbola. With a little algebra 
that's too long to show here, we can 
prove that the hyperbola opens along 
the y-axis, is centered at the point (1/2, 
1/2), and has limiting values of p given 
by: 

p = \{\±i^2) (35) 



the equation for the hyperbola because 
we don't really need it. We can simply 
solve Equation 34 for y, just as we did 
for Equation 29. The equation has the 
form of a quadratic with: 

a = -1; b = 1; c = 2p (p - 1) (36) 
The quadratic formula gives us: 



Using Equation 27 as before, this I'm not showing the precise form of 



K = !(l±Vl-80(0-l)) (37) 

This curve is plotted as the heavier line 
in Figure 5. Between the two curves, 
defined by u = and v = 0, we've 
defined the boundaries of all four 
cases. You might be wondering what 
happened to Case 3, since it's not 
shown on the curve. As it turns out, 
this case is represented by a tiny sliver 
where the two curves intersect near the 
point p = 0.85, y = 0.5. The region is 
too small to show up on this graph. Its 
precise location and shape turn out to 




Here? 



Here? 
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FIGURE 5 

Boundaries. 




be immaterial anyway because, while 
we've shown the graph for positive 
and negative values of B and y, the 
restrictions given by Equations 1 3 and 



14 demand that fi and y be negative. 
Thus, we're really only concerned with 
the lower left-hand quadrant of 
Figure 5. 



WHICH REGION? 

Only Cases 1 and 4 satisfy the inequal- 
ity in Equation 24. Within that lower 
left quadrant of Figure 5, we are 
allowed to operate in the triangular- 
shaped part of Case 1 and the part of 
Case 4 curving in from the left. 

I'm going to assert now that only the 
Case 1 portion will do us any good. 
Here's why: Consider what happens as 
we start getting closer to convergence. 
Presumably, our boundary points Pj 
through P 3 will get closer and closer 
together, and the segment of the func- 
tion they define will become more 
nearly a straight line. When the curve 
is a straight line, the quadratic of 
Equation 5 degenerates as the quadrat- 
ic coefficient C approaches zero. 
Let's see what must happen to B and y 
as this degeneration takes place. 
Setting C to zero in Equation 1 1 gives 
us: 

-y, + 2y 2 - y 3 = 



}ou ever think being an engineer could be so 
" You've worked hard to acquire the 
experience and expertise you need to push the 
i tecruiological envelope. 

- Yet instead of solving your most challenging 
^pblerrjs, |b)i [ can spend most of your day bump- 
ing into dead ends. 

Ipfhe best use ofyourtime, oryourmind. 
> why HP is offering 



r way for you and 
ur team to work. Ourphi- 
ophy is simple: give you 
"hforniation you 
ieed, in a relevant and use- 
ale format, so you can find 
oblems through logical 
Wdng;not guesswork. 

How HP's design philosophy gives you 



To spend 
less time guessing, 
start here. 



ter insight. 

.y^need to', have the right tools available, 
e ve deyejoped aicomplete range of affordable 
•ftware tools - from ROM monitors to Background 
" leandru^perfoimanceemulators-sono 
tt problem you're, troubleshooting, apply- 
t solution is easy. 




Second, you need to have information that's not 
onlyinalanguage you understand, butinacontext 
that's relevant. That's why real-time, high-level 
language debugging is our standard. Because if 
you can see how your code relates to the system - 
right when an error occurred - you'll immediately 
know what caused the problem. 

Finally, you need tools that will get used. We 
designed easy-to-use,open 
systems with verified connec- 
tions to leading software 
vendors. So you'llfeel confi- 
dentchoosingthebesttools, 
knowing that if they work 
well together, so can you 
and your project team. 

Get started today. 

For faster insight into your software design prob- 
lems, call 1-800-447-3287, Ext. 8369, and ask 
for our free "Your Solution's In Sight" Kit. It in- 
cludes a product demonstration disk, a Software 
Designer Concept brochure, and technical litera- 
ture on HP's entire family 
of software solutions. 

Or just turn the page. 
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or the following: 
2p-y = l 



(38) 



This, of course, is a straight line, which 
I've plotted as the dashed line in Figure 
5. As we are iterating and getting clos- 
er to the true root, we can expect the 
values for each iteration to migrate 
toward this line; the actual point will 
depend on the slope of the function 
curve. In the lower-left quadrant, this 
line lies wholly in the Case 1 region. 
Thus our point can indeed migrate to 
the straight line. The Case 4 region, on 
the other hand, is cut off from the solu- 
tion line by the Case 2 region sur- 
rounding it. Because of this, starting a 
quadratic step from within Case 4 is a 
fruitless exercise; we'll end up having 
to switch back to bisection anyway. 
Case 4, in fact, corresponds to the case 
of an acceptable, but poor, fit similar to 
that of Figure 4. 

What we've gained out of all these 
equations and graphs, then, is the con- 
viction that only the small triangular 
region in the lower right corner of the 
lower left quadrant will be suitable for 
inverse quadratic interpolation. In all 
other regions, we should stick with 
bisection. 

Armed with this conviction, we're 
finally prepared to define a useful cri- 
terion for accepting the interpolated 
root. Look again at the four cases; the 
only difference between Case 1 and 
Case 3 is the sign of u. If we had to 
worry about falling into Case 3, we 
would have to test the signs of both u 
and v. Since Case 3 only lies in that 
tiny sliver in the upper right quadrant, 
we needn't worry about it. We can 
never fall into that region. We can 
ignore the sign of u and accept as our 
criterion the condition: 



v>0 



(39) 



The value of v is given by Equation 
22. This leads us to a simple test for 
acceptable interpolation: 

>'3(y3-v 1 )-2v 2 (y 2 -y,)>0 (40) 
Having struggled with all those equa- 



tions and graphs, you can now forget 
them. They served only to establish the 
condition in Equation 40. 

THE ALGORITHM 

Now we can give the complete algo- 
rithm for Method X: 

• Given two points x, and x 3 straddling 
a root, bisect the region: 



x 2 



X\ + x 3 
= 2 



(41) 



and evaluate the function to get y 2 . 

• Test the sign of y 2 y3- If it is negative, 
exchange the points P, and P 3 to make 
it positive. 

• Test the sign of: 

v= y 3 (y 3 -y,)-2y 2 (y 2 - V| ) 

If the sign is negative, replace P 3 by P 2 
and repeat. If the sign is positive, com- 
pute the coefficients: 



B = 



x 2 ~ x \ 



c _ -y\+ 2y 2 - y 3 
(y 3 - yaXys - y,) 



• Compute the estimated root from: 
x = *j - #yi(l - Cy 2 ) 

and evaluate the new function value, 

y = f(x). 

• Test the sign of y to determine which 
points to use for the next cycle. If yy, 
< 0, replace P 3 by (x,y). Otherwise 
replace P) by (x,y) and P 3 by P 2 . 

• Repeat until two successive roots 
are equal within a specified error, 
epsilon. 

IMPLEMENTATION 

Now that we finally have the algo- 
rithm, implementing it in code is 
straightforward enough. The results are 
shown in Listing 1. I tested the algo- 
rithm using our now-standard test 
function, which is shown on the fol- 
lowing page: 



LISTING 1 

Second-order iteration. 

enum {none, bad.data, no.convergence} 
error; 

// iterative root-finder using 

// alternating bisection and inverse 

// quadratic interpolation 

// call as: x2 = iterate(xl,x3,f, 
// eps.imax.error); 
// where xl,x3: initial points 
// bracketing root 

// f(x): name of user-supplied function 

// eps: accuracy required in x 

// imax: maximum number of iterations 

// function return value contains result 

// error: output flag 

// error=none -> normal return 

// error=bad.data -> initial points 

// don't bracket root 

// error=no.convergence -> cannot 

// converge in inax tries 

double iterate(double xl, double x3, 
double (*f) (double) .double eps.int 
imax.int ft error){ 
double x2,yl,y2,y3,b,c ) temp,y21,y31, 
y32,xm,ym; 
error=none; 
yl=f(xl); 

if (yl=0.0) return xl; 
y3=f(x3); 

if(y3=0.0)return x3; 
if(y3*ylX>.0){ 

error=bad_data; 

return xl; 

} 

for(int i=0;i<imax;i++){ 
x2=0.5*(x3+xl); 
y2=f(x2); 

if(y2=0.0)return x2; 
if(abs(x2-xl)<eps)return x2; 
if(y2*ylX>.0){ 

temp=xl; 

xl=x3; 

x3=temp; 

temp=yl; 

yi=y3; 
y3=temp; 

y21=y2-yl; 
y32=y3-y2; 
y31«y3-yi; 

if (y3*y31<2. 0*y2*y2lH 



x3=x2; 
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y3=y2; 

else{ 
b=(x2-xl)/y21; 
C=(y21-y32)/(y32*y31); 
xm=xl-b*yl*(1.0-c*y2); 
ym=f (xm) ; 

if(ym=0.0)return xm; 
if(abs(xm-xl)<eps)return xm; 
if(ym*yl<0.0){ 

x3=xm; 

y3=ym; 

} 

else{ 
xl=xm; 
yl=ym; 
x3=x2; 
y3=y2; 

} 

} 



TABLE 1 

Iterative results. 



} 



error=no_convergence; 
return x2; 



Iterations 


Estimated Root 


Error 


1 


4.917669 


1.515485 


2 


3.214061 


0.188123 


3 


3.398273 


0.003911 


4 


3.402173 


0.000011 




3 402184 






fix) = x-4sinx + e- xl6 -5 

The results are shown in Table 1 . As 
you can see, the speed with which the 
algorithm locks in on the root, once it's 
gotten in the neighborhood, is impres- 
sive indeed. William Press writes that 
Brent's method is the method of choice 
for general-purpose use. Although 
either Brent's method or Method X 
should give equally good performance, 
since they use the same philosophies, I 
claim that Method X is the method of 
choice because the equations are a tad 
simpler. Use this function in good 



health. You won't find a better one any- 
where in the known universe. I 



Jack W. Crenshaw is a staff scientist at 
In Vivo Research in Orlando, FL. He 
has developed numerous analysis and 
real-time programs and holds a PhD 
in physics from Auburn University. 
Crenshaw can be reached via e-mail at 
72325. l327@compuserve.com. 
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Looking at information one piece at a 
time can send you down a lot of blind 
alleys. That got us thinking. What if an 
emulator could help you go straight to 
the answer? 

With dual-ported emulation memory and 
a choice of foreground or background 
monitors, you could make measurements 
in real-time - without interrupting your 
design process. 

If you could do high-level C/C++ dynamic 
debug and analysis, you'd spend less time 
deciphering code. 

Common debugger interfaces would 
make sharing data with your team easier. 
And if you could get all this on an em- 
ulator for virtually any processor you 
have, your problems would be solved. 

If that sounds like a faster approach to 
you, call 1-800-447-3287, Ext. 8369 
for a free demonstration disk Or ask for 
an HP engineer. 

And speak to us directly. 

There is a better way. 
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